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^ ABSTRACT 
(N 
l> 

' We model the photoelectric emission from and charging of interstellar dust and ob- 

Q> . tain photoelectric gas heating efficiencies as a function of grain size and the relevant 

c~| . ambient conditions. Using realistic grain size distributions, we evaluate the net gas 

heating rate for various interstellar environments, and find less heating for dense re- 
O I gions characterized by Ry = 5.5 than for diffuse regions with Ry = 3.1. We provide 

I fitting functions which reproduce our numerical results for photoelectric heating and re- 

' combination cooling for a wide range of interstellar conditions. In a separate paper we 

will examine the implications of these results for the thermal structure of the interstellar 
^ ■ medium. Finally, we investigate the potential importance of photoelectric heating in 

^ . H II regions, including the warm ionized medium. We find that photoelectric heating 

could be comparable to or exceed heating due to photoionization of H for high ratios 
of the radiation intensity to the gas density. We also find that photoelectric heating 
by dust can account for the observed variation of temperature with distance from the 
galactic midplane in the warm ionized medium. 

Subject headings: dust — extinction — HII regions — ISM: clouds 



1. Introduction 

Photoelectric emission from dust grains dominates the heating of diffuse interstellar gas clouds 
as well as the photodissociation regions at the surfaces of molecular clouds. This mechanism 
therefore plays an important role in the dynamical evolution of the interstellar medium, and in 
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shaping the hne emission spectra of galaxies. Heating by photoelectric emission from interstellar 

grains was first considered by Spitzer (1948), and there have since been a number of reassessments 
of the heating rate (Watson 1972; deJong 1977; Draine 1978; Tielens k Hollenbach 1985; Bakes & 
Tielens 1994, Dwek k Smith 1996). 

For a given grain, the heating rate depends on the grain size, composition, and charge state, 
as well as on the spectrum of the illuminating radiation. Because photoelectric yields are enhanced 
for small grains (Watson 1972, 1973), estimates of the net photoelectric heating rate in interstel- 
lar gas are sensitive to the adopted grain size distribution, which should be consistent with the 
observed extinction curve (which shows strong regional variations) as well as with the observed 
emission from interstellar dust grains, from the near-infrared to the microwave. Previous estimates 
for photoelectric heating rates did not always use grain size distributions which were consistent 
with these constraints. While the measured extinction curve by itself does not suffice to uniquely 
specify the grain size distribution, in the present study we will consider size distributions which are 
consistent with the observed extinction in different regions, with either the minimum or maximum 
permissible population of ultrasmall grains (Weingartner & Draine 2000). 

In addition to using size distributions consistent with observations, we also model the photo- 
electric emission process and associated grain charging in detail, using realistic yields for graphitic 
and silicate grains, a realistic distribution of photoelectron kinetic energies, and new estimates for 
electron sticking efficiencies for small grains. The resulting photoelectric heating rates are calcu- 
lated for grain size distributions consistent with extinction curves characteristic of diffuse clouds 
{Ry = Av/Eb^v ~ 3.1), intermediate density regions {Ry ~ 4.0), and dense clouds {Ry ~ 5.5). 

In §2 we discuss our adopted model for photoelectric emission from grains, expanding on the 
treatments of Draine (1978) and Bakes k Tielens (1994). In §3 we use recent experimental results 
to obtain improved estimates for the electron sticking coefficient, as a function of grain size; our new 
sticking coefficients differ significantly from previous estimates. In §4 we characterize the relevant 
ambient conditions (which are important in grain charging) and in §5 we evaluate the photoelectric 
heating efficiency as a function of grain size. We obtain net heating rates in H I regions in §6 and 
in §7 we investigate the potential contribution of photoelectric heating in H II regions, finding that 
grain photoelectric heating may explain observed high temperatures in the warm ionized medium. 
We briefly summarize our results in §8. 



2. Photoelectric Emission from Grains 

2.1. Grain Characterization 

We consider both graphitic and silicate grains. We adopt the description of graphitic grains 
given by Li k Draine (2000), in which the smallest grains are polycyclic aromatic hydrocarbon 
(PAH) molecules, the largest grains consist of graphite, and grains of intermediate size have optical 
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properties intermediate between those of PAHs and graphite. We characterize the graphitic grains 
by an effective radius a, the radius of a sphere containing the same number of C atoms, in which C 
atoms contribute a mass density p = 2.24 gcm^'^, the density of ideal graphite. The smallest grains 
certainly are not spherical and probably are not well characterized by the ideal graphite density; 
in this case the effective radius simply indicates the number Nc of C atoms in the grain: 



In following sections, the Nq-q relation (1) will be used to incorporate experimental data on 
hydrocarbon molecules into our description in terms of the effective radius. Since small molecules 
that extrapolate to silicates as a — ^ oo have not been studied, we will not need a similar relation 
for silicates. In this case, a is simply the radius of a spherical silicate grain. 



The electrons in a solid are confined in a potential well (see, e.g., Lang & Kohn 1971). For 
neutral bulk material, the work function W is the energy of a free electron with zero kinetic energy 
minus the energy of an electron in the highest occupied state in the solid. (Throughout this 
section, refer to Figure 1 for schematic depictions of the electron-confining potential.) In this case, 
the ionization potential IP, i.e., the difference in energy between infinity and the highest occupied 
state, equals W. 

For small grains, IP ^ W. We first consider grains which are not so small that quantum 
effects shift the energy levels from their bulk values. If Z > 0, then an electron that is liberated 
from the grain material feels a Coulomb attraction to the grain, which now has charge {Z + l)e. 
For a conducting sphere, IP — W can be determined by calculating the work needed to remove 
the electron from a point just outside the sphere to infinity (including the term due to the image 
charge as well as that due to the net charge on the sphere) and subtracting the corresponding work 
associated with a plane surface. The result is that IP = W + {Z + l/2)e^/a (see, e.g. Makov, 
Nitzan, &; Brus 1988; Seidl &; Brack 1996).^ For nonconducting materials, the term e^/2a above 
should be multiplied by (e — l)/e, where e is the dielectric function (Makov et al. 1988). This 
correction is fairly small for candidate grain materials, so we will ignore it. 

When Z < the "extra" electrons, which occupy the "LUMO" (for "lowest unoccupied 
molecular orbital") of the neutral, are separated from the top of the valence band by the band 



^Surprisingly, there has been some controversy in the hterature as to the value of the constant added to Z. See 
Makov et al. (1988), Moskovits (1991), Seidl & Brack (1996). 
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2.2. Electron Energy Levels 
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gap, with extra energy Ey^g (assuming the valence band is full in the neutral). The electron affinity 
EA{Z) is the difference in energy between infinity and the LUMO for the grain of charge [Z — l)e 
created by adding the electron. Thus, IP{Z < 0) = EA{Z + 1). We also define the "valence band 
ionization potential" IPv{Z), equal to the difference in energy between infinity and the top of the 
valence band. Of course, IPv{Z) = IP{Z) when Z > 0. 

Since it is not yet possible to calculate the quantum shifts in energy levels for very small 
grains, we adopt an empirical approach. Figure 2 shows the first and second ionization potentials 
/P(0) and -f-P(l) for various aromatic hydrocarbons, including CeHg benzene, C24H12 coronene, 
and C38H20 benz[42]. The polycyclic aromatic hydrocarbons (PAHs) of interest here appear to be 
fit by 

,^(z).vr+(z + l)^ + (z + 4(^iA) , (2) 

with the graphite work function W = 4.4 eV . This estimate for IP{Z) is plotted in Figure 2 
for Z = and 1. For the first ionization potential IP{0), equation (2) is numerically very close 
to the estimates of Bakes (1992) and Bakes and Tielens (1994), but for IP{1) our estimate is 
~1 eV below their estimate for PAHs with ~15-30 C atoms. Note that equation (2) reproduces 
the classical electrostatic result for large a. Since experimental evidence is not available for other 
ionization states, wc will adopt equation (2) for all values of Z. 

For silicates, few laboratory experiments are available to guide us. We adopt W = 8cV , since 
the photoelectric yield for lunar dust drops steeply around that energy (Feuerbacher et al. 1972). 
Nayak et al. (1998) have made theoretical calculations of the ionization potentials of small silica 
clusters (Si02)n, with n = 1 to 6. They find IP{0) = 11.19 eV (16.08 eV) for n = 2 (5), and 
intermediate results for other values of n. If we employ equation (2) with = 8eV and a ^ 3.2 A, 
then we find IP{0) ~ 12 eV , which is in the Nayak et al. range for small clusters. Thus, we adopt 
equation (2) for both carbonaceous and silicate grains, with W = 4.4 and 8 eV , respectively. 

We expect that 

EA{Z <0) = W-E^g{a)+(^Z-^^^ + O{a-^) , (3) 

where the presence of an 0(a^^) term is expected since such a correction was required to fit the 
experimental results for IPv{Z > 0) (eq. 2). The "infinite limit" for PAHs - graphite - is a 
"semimetal" , with a slight overlap between the top of the valence band and the conduction band 
(Spain 1973). For metals and semimetals, ii'bg = 0. Many silicates begin to absorb strongly around 
2500-2000 A, corresponding to a band gap of 5-6 eV ; we will adopt £^bg = 5eV, independent of a, 
for silicates. 

Values of EA{0) are shown in Figure 3 for selected neutral PAHs, benzene, and Ceo- The 
relation 

(._^)f!_f!(_lA_) „ 
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is also plotted in Figure 3, for Z = and W = 4.4 eV. For a — > oo this has the theoretically 

expected behavior, and it provides a reasonable fit to the experimental values of EA{0) for benzene 
and the PAHs in Figure 3 (for which a « 2.3 — 4. 3 A). We adopt this equation (4) for carbonaceous 
grains. For silicates, we adopt 

EA{Z) = W-E^^+(^Z^^y^, (5) 

with W — E\^g = 3eV . This yields EA{0) « 0.73 cV for a f» 3.2 A, which is somewhat low compared 
with the values of 2-3 eV inferred for small silica clusters (Wang et al. 1997). 



2.3. Photoelectric Yields 

The photoelectric yield is the probability that an electron will be ejected following the absorp- 
tion of a photon. Calculation of the yield from first principles is not yet possible. Draine (1978) 
displayed yield measurements from experiments on bulk samples of materials of interest, includ- 
ing graphite, lunar dust, and silicon carbide. Watson (1972, 1973) pointed out that the yields of 
submicron particles are expected to be enhanced relative to the bulk yields, because of the finite 
electron escape length. An electron excited somewhere in the vohime of the sample can lose energy 
during its journey to the surface, through interactions with electrons and with phonons. The result 
is that the fraction of electrons which escape energy loss goes roughly as exp(— x/^e), where x is 
the distance the electron has traveled and Ig is the "electron escape length." Generally, the photon 
attenuation length, (equal to the e-folding length for the decrease in radiation intensity as it 
propagates into the material), exceeds le, and the electrons excited deep inside a bulk sample do 
not reach the surface. Small grain sizes limit the distance from electron excitation to surface, so 
the yield is enhanced. 

A detailed and accurate model for the above geometrical effect is not yet available. A series of 
experiments (Schmidt-Ott et al. 1980; Burtscher & Schmidt-Ott 1982; Burtscher et al. 1984; Miiller 
et al. 1988) on free silver spheres with radii between 27 and 54 A found yield enhancements in excess 
of expectations by factors of several. A later theoretical effort (Faraci, Pennisi, & Margaritondo 
1989) was able to reproduce the results, but only by assuming that the condition for an electron 
to escape, when incident on the surface from within, depends only on its energy and not on its 
direction of motion. The estimates adopted here for size-dependent yields from grains are based on 
Watson's (1973) simplified model for the geometrical effect and must be regarded as provisional. 
Reliable calculations of effects involving photoelectric emission from interstellar dust will not be 
possible until experiments have been performed on submicron grains of appropriate composition. 

For negatively charged grains we will distinguish between "photoelectric" ejection of electrons 
from the valence band and "photodetachment" of excess "attached" electrons. 
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2.3.1. Valence Band Electrons 



We assume that, when Z > 0, the highest occupied energy level is very close to the top 
of the valence band, since the number of electrons which have been removed from the grain is 
small compared with the total number of electrons in the grain. Thus, for Z > 0, the threshold 
photon energy for photoelectric emission is given by hfpetiZ) = IPy{Z) = IP{Z). When Z < — 1, 
hiypct{Z) > IP{Z), because the electron has to overcome the repulsive Coulomb barrier. Suppose 
the tunneling probability becomes significant when the electron energy exceeds the potential at 
infinity by .Emm- Then hu^^^^Z) = IPv{Z) + £^mm- Thus, we take 



{IPv{Z,a) 



{ IPv {Z,a) + E^i^ {Z, a 



Z > 
Z < 



(6) 



We estimate Eiam{Z,a) as the energy for which the tunneling probability (evaluated using 
the WKB approximation) equals 10~^; the tunneling probability rapidly increases as the excited 
electron energy increases above this value. The barrier consists of the Coulomb potential —{Z + 

l)e^/r and the image potential — e^a^/2r^(r^ — a^), where r is the distance from the center of the 
spherical grain. ^ We adopt the following expression for .Emin, which reproduces the results of the 
WKB calculation fairly well when -(Z + l)e^/a > 0.5 eV : 



E^iniZ <0,a) = -{Z + l)- 



l+ 



27 A 



0.75' 



(7) 



When —{Z + l)e^/a < leV, equation (7) overestimates E^, 
severity as —{Z + l)e^/a — 0. However, in this limit Emm - 
will not substantially affect the photoemission calculations. 



n as defined above, with increasing 
as well, so that this overestimate 



We adopt a very simple model for estimating the photoelectric yield Y{hi', Z,a). We first 
consider the yield t/q of electrons that traverse the surface layer, emerge from the grain surface, and 
have enough energy to overcome the image potential (we will call these "attempting" electrons). 
For Z < 0, every attempting electron will escape, but when Z > 0, some of these electrons have 
insufficient energy to escape to infinity and instead fall back to the grain. We assume the following 
form for yg: 

yo = yo{Q)yi{hiy,a) , (8) 



where the parameter is given by 



e = 



hu 



hup^t + {Z + l)eVa 



Z>0 
Z <0 



(9) 



^Within a few A of the surface, the image potential deviates from the classical expression, and "saturates" to a 
constant value (see, e.g., Lang & Kohn 1971). However, this saturated portion of the potential is always classically 
allowed, and therefore does not affect the WKB estimate of the tunneling probability. 
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and the factor yi{hu,a) accounts for the size-dependent geometrical yield enhancement discussed 
above; yi depends on hv because the photon attenuation length does. 

We assume a parabolic energy distribution for the attempting electrons:^ 

f^{E) = 6(^~^'ow)(-^high--g) ^^Q^ 
(-E^high — -E^low)^ 

where f%{E)dE gives the fraction of attempting electrons with energy (with respect to infinity) 
between E and E + dE. When Z < 0, E'lo^ = -Emin (eq- 7) and -Ehigh = -E'min + hu — /if pet; when 
Z >0, E'low = —{Z + l)e^/a and E'high = hv — hv^et- 

Let 2/2 be the fraction of attempting electrons which escape to infinity: 

y2{hv, Z, o) = I ■^0^''""' ^^"^E^-^) = ^high(^high - 3£;iow)/ (^high - -Eiow)^ , ^ > ^^^^ 

Our resulting expression for the photoelectric yield is 

Y{hv,Z,a) =y2{hv,Z,a)mm.[yQ{Q)yi{a,hv),l] . (12) 

Draine (1978) found that the following estimate for yi reproduces Watson's (1973) detailed 
results based on Mie theory, to within 20%: 

_ //3y a^-2Q + 2-2exp(-a) 

\aj /?2 - 2/3 + 2 - 2exp(-/3) ' ^ ^ 

where /? = a/la ^ = o/^a + «/^e- 

Electron escape lengths for most materials reach a minimum in the vicinity of tens of eV and 
increase at higher and lower energies. The data are meager for electron energies below 10 eV , but 
Martin et al. (1987) found, for thin carbon films, that Ze fell into a broad minimum of 6 A around 

40 cV and rose to 9 A at 6eV. McFeely et al. (1990) found, for Si02, that le varied from 5.7 to 
6.8 A for electron energies between 8 and 20 eV , and one might expect similar values for silicates. 
Bakes &: Tielens (1994) assumed 1^ = 10 A, independent of energy, and we use the same value here, 
for both graphite and silicates. 

The photon attenuation length is given by 

la= , , , (14) 

47rlm(m) 

where A is the wavelength in vacuo and ra(X) is the complex refractive index. Graphite is a highly 
anisotropic material, so that the dielectric function is a tensor. This tensor may be diagonalized 



parabolic form is consistent with typical laboratory photoelectron distributions which drop to zero when E = 
and E = hv — W and peak somewhere in between (Draine 1978). 
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by choosing Cartesian coordinates with two of the axes lying in the "basal" plane; the third axis, 
normal to the basal plane, is called the "c-axis" . For graphite, we take 



7-1 - ^ 

« " A 



2 1 
-Im(m_L) H — Im(m|| 

3 3 



(15) 



where m± and my are for the electric field perpendicular and parallel to the c-axis, respectively. 
In computing we use dielectric functions from Draine & Lee (1984) and Laor & Draine (1993), 
modified to remove a silicate feature of crystalline origin that is not present in the observed inter- 
stellar extinction or polarization (see Wcingartncr & Draine 2000). Since PAH dielectric functions 
are not available, we use graphite dielectric functions for all carbonaceous grains. 

The final step in our yield prescription is to estimate ?/o(0) for carbonaceous and silicate grains. 
For carbonaceous grains, we follow Bakes &: Tielens (1994) and choose yo such that Y^hv, Z = 0,a = 
3.7 A) approximately reproduces the photoionization yield for coronene, as measured by Verstraete 
et al. (1990). This gives 

,„(e)= 9xio-3(e/HO^ 

* 1 + 3.7 X 10-2 (e/VF)^ ^ ' 

with W = 4.4 eV (see Figure 4). Draine (1978) shows, in his Figure 1, that the bulk graphite yield 
measured by Fcucrbachcr &: Fitton (1972) is unusually low, and lies more than an order of magnitude 
below the yield for anthracene (Fujihira, Hirooka, & Inokuchi 1973). It is important to note that 
our equation (16) gives yields that substantially exceed those measured for bulk graphite.^ This 
apparent inconsistency could imply one or more of the following: 1. The yield for bulk graphite 
differs substantially from that for bulk coronene. 2. The measured yields for bulk graphite are 
wrong. 3. Our prescription for yi^a^hf) is wrong. Clearly, the adopted yields for carbonaceous 
grains are highly uncertain. 

The situation for silicates is no better. Feuerbacher et al. (1972) measured the yield for a 
powdered sample of lunar dust, and found that Y decreases rapidly as photon energy is decreased 
below 14 eV . However, they noted that the use of the powdered sample might result in yields that 
are too low. Thus, we choose yo for silicates such that the bulk yield somewhat exceeds the results 
of Feuerbacher et al. at 14 eV , but does not drop as rapidly for lower values of hv: 

0.5 (G/M/) 

^"(Q)= l + 5(9/^ ' ^^'^ 

with W = Se\ . 

In Figure 5 we plot Y{hv) for neutral carbonaceous and silicate grains of several sizes. 



''For example, at hv = lOeV , where Feuerbacher & Fitton (1972) measured a yield of 8 x 10~^ for bulk graphite, 
our equation (16) gives yo = 2.7 x 10~^, a factor 3.4 larger than the laboratory value. For comparison, the Bakes & 
Tielens (1994) expression gives yo = 7.8 x 10~^, an order of magnitude larger than the experimental value. 



-9- 



2.3.2. Photoelectric Ejection of Attached Electrons (Photodetachment) 

When Z < 0, the —Z attached electrons occupy energy levels above the valence band, if the 
latter is full in the neutral. We take the photodetachment threshold energy /ifpdt to be 

hupAtiZ < 0) = EA{Z + 1, a) + E^i^iZ, a) (18) 

with electron affinities EA given by equations (4) and (5), and E^nin by (7). Bakes & Tielens 
(1994) neglected £^min(Z, a) and took /it'pdt = EA{Z + l,a). We assume that E = hv — hv^^i for 
photodetached electrons, since attached electrons lie in a narrow range of energies. 

We assume that an oscillator strength /pdt is associated with photodetachment transitions 
to the continuum. The photodetachment cross section f^pdt(^^) has been measured for CgFg 
(Christophorou, Datskos, & Faidas 1994). The measured cross section has considerable structure, 
but can be roughly approximated by 

. ^ 27re^fe/pdt X 
aMZ,a) = -Z 3^^^^^ (1 + ^2/3)2 ' (19) 

where x = {hu — /iz/pdt)/A£^ and the peak in (jpdt occurs at hv = /ifpdt + We take ^E = 3eV 
and oscillator strength /pdt = 0.5; thus 

apAi{hv,Z,a) = 1.2 X \{)-^^ cvc?\Z\-^^-^^^j^ for Z < 0. (20) 



2.4. Grain Charging 

Since the photoemission depends on the grain charge, it is necessary to know the distribution 
of charge states for the grains. In statistical equilibrium, 

/z(Z)[Jpe(Z) + Jio„(Z)]=/z(Z + l)Je(Z+l) , (21) 

where fz{Z) is the probability for the grain charge to be Ze, Jpe is the photoemission rate, Jion is 
the positive ion accretion rate, and Jg is the electron accretion rate. 

Bakes k, Tielens (1994) discuss the most positive and most negative charges that a grain could 

possibly acquire (Zmax^ and ZminC, respectively). The most positive charge is one proton charge 
more than the highest charge for which an electron can be ejected, i.e. for which hvp^t = IP < 
hvmax, the maximum photon energy in the radiation field (= 13.6 eV , in an H I region). Thus, 



■^ma.x — int 



fet^max -w a 1 _ a3A\ / o.sA 

14.4eV A 2 a " 



(22) 



where int [x] denotes the greatest integer less than x. The minimum allowed charge Z^me is the 
most negative charge for which autoionization does not occur. Bakes & Tielens (1994) assume 
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that this requires EA{Zea + 1) > 0. However, large grains can be charged more negatively 
without emitting a substantial electron current when the tunneling probability is low. We take the 
autoionization threshold potential Uait to be that at which the electron current is ~ lO^^s^"*^; we 
evaluate the tunneling probability using the WKB approximation and assume an attempt frequency 
of f=i 2 X 10«s-i(cm/a). We find 

-Unit _ / 3.9 + 0.12 (a/ A) + 2 ( A/a) , for carbonaceous 
V ~ 1 2.5 + 0.07 (a/ A) + 8 ( A/a) , for silicate ' ^ ' 

this agrees very well with Draine & Salpeter's (1979) simple estimate of the potential at which field 
emission becomes important. The most negative allowed charge is then given by 



Zm\n — int 



14.4V A 



+ 1 . (24) 



For carbonaceous grains, equations (23)-(24) give Zmin = for a < 2.92 A, Z^ai-a = —1 for 2.92 < 
a < 5.75 A, Zmin = -2 for 5.75 < a < 8.40 A. For silicate grains, Z^i^ = for a < 2.40 A, 
^min = -1 for 2.40 < a < 6.96 A, Zmin = -2 for 6.96 < a < 10.8 A. By iteratively applying 
equation (21) and normalizing, fz can be found for all Z. 

The photoemission rate is given by 

Jpe = 7ra^ / dvYQ^y^sW)—^ + / dva^AtW)-r^ , (25) 

where Qabs is the absorption efficiency, Ui, is the radiation energy density per frequency interval, 
and c is the speed of light. The second term in equation (25) is only present when Z Bakes & 
Tielens (1994) assumed Qabs a, valid for grains with a < 100 A. Since we consider larger grains, 
we cannot adopt this approximation. We evaluate Qabs for graphitic grains using the prescription 
of Li & Draine (2000)^ and for silicate grains using a Mie theory code derived from BHMIE (Bohren 
&; Huffman 1983) with dielectric functions given by Draine Sz Lee (1984) and Laor Sz Draine (1993), 
but modified in the ultraviolet (Weingartner &; Draine 2000). 



The accretion rates are given by 



1/2 



MZ)=nMZ){ 7raV(ri,ei) , (26) 

where rij is the number density of species i, the sticking coefficient Si{Z) is the probability that 
species i will transfer its charge if it reaches the surface of a grain of charge Ze, rrii is the particle 
mass, T is the gas temperature, and J is a function of = akT/qi^ and = Ze/qi {qi is the charge 
of species i and k is the Boltzmann constant). Expressions for J can be found in Draine & Sutin 
(1987). The electron and ion sticking coefficients Sg and Sj are discussed below in §3. 



^Li & Draine (2000) give optical properties for neutral and ionized grains; we use the former when Z = Q and the 
latter when Z ^ 0. The distinction is not critical, as in the vacuum ultraviolet there is little difference between the 
ionized and neutral grain absorption cross sections. 
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3. Sticking Coefficients 
3.1. Electrons 

A low-energy electron impinging on a macroscopic solid surface has some probability pes of 
elastic scattering. We will assume that pes ~ 0.5, so that the maximum possible value of the 
sticking coefficient Sg would be (1 — Pes) ~ 0.5. 

3.1.1. Attachment to Neutral Grains 

We first consider electron attachment to neutral grains. If the electron affinity EA > [as 
it is for a > 3.5 A, which Weingartner & Draine (2000) take to be the lower cutoff for the grain 
size distribution], then the approaching electron accelerates due to its polarization of the grain, 
arriving at the surface with a kinetic energy of order EA. Even if it enters the grain material [with 
probability (l—pes)], the electron may fail to undergo an inelastic scattering event, in which case it 
passes through the grain and returns to infinity. The probability of undergoing inelastic scattering 
is approximately (l — e^'*/''=), where /g ~ 10 A, the "electron escape length" (discussed above), is 
roughly the mean free path against inelastic scattering within the grain material. 

If the impinging electron docs undergo inelastic scattering, then it transfers some of its energy 
to internal degrees of freedom of the grain; if the energy so transferred exceeds the initial kinetic 
energy (~ 2kT), then the electron is trapped by the grain potential. However, until the grain 
radiates this energy away, there is a nonzero probability per unit time that the internal degrees 
of freedom of the grain will transfer enough energy back to the electron to eject it. Let Prad be 
the probability of "radiative stabilization" of the negatively charged grain, i.e. the probability of 
radiating away ~ kT of energy before the electron is ejected. We expect prad 1 for macroscopic 
grains, but prad < 1 foi' grains with a small number of internal degrees of freedom. The electron 
attachment sticking coefficient can then be written as the product of three factors: 

Se(Z = 0)«(l-pe.)(l-e""/^^)prad • (27) 

The stabilization probability Prad can be estimated from the electron affinity EA, the density 
of states of the excited negatively charged grain created by the electron attachment, and the 
probability per time of photon emission from the excited grain (Allamandola, Tielens, & Barker 
1989; Tielens 1993). However, since these parameters are poorly known, we will instead adopt an 
empirical approach. 

In Figure 6 we show sticking coefficients Se = {av) / [{8kT / nmey^^Tra'^ J] for various hydrocar- 
bon molecules, where J{akT /e^ ,0) is given by Draine k. Sutin (1987) and the rate coefficients {av) 
for electron attachment are determined experimentally. Also shown in Figure 6 is a semiempirical 
fit: 

^e(^ = 0) = 0-5(l-^"°^'-) l + e(^o-iVc) ' (28) 
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where A^c is the number of atoms other than H in the molecule, and Zg = 10 cm. 

The measured sticking coefficients for Nq < 20 are in reasonable agreement with our semiem- 
pirical fit (eq. 28). It is notable that Tobita et al. (1992) found that in many cases the sticking 

coefficient increases when the electron energy increases. For example, the sticking coefficients for 
tetracene C18H12 and perylene C20H12 would be larger by an order of magnitude for gas temperature 
T ~ 2000K, bringing them up to or even above the values given by equation (28). 

The measured Sg for Cgo and C70 at T < 500K are remarkably low, probably due to the 
unusual rigidity and symmetry of these molecules. For Cqq, the symmetry forbids capture of 
s— wave electrons, and the centrifugal barrier to p— wave electrons results in an apparent activation 
energy of 0.26eV (Tossati & Manini 1994). At high temperatures the sticking coefficients for Ceo 
and C70 are close to the values given by equation (28). 

We adopt equation (28) for Se{Z = 0) for silicates as well as for carbonaceous grains, since no 
laboratory data are available for silicates. In this case, Nc is just a parameter defined by equation 

(1)- 

3.1.2. Attachment to Negatively Charged Grains 

As a grain acquires more electrons, the electron affinity decreases. We take the sticking coeffi- 
cient to be zero when Z = Zj^in, since the more negatively-charged state, even if it formed, would 
autoionize. Thus, 

SeiZ<0) = \'/^ = ^^ 'l^l^"" . (29) 

3.1.3. Recombination with Positively Charged Grains 

Electrons arriving at positively charged grains are expected to recombine provided that (1) 
they do not reflect elastically from the surface and (2) they are able to scatter inelastically before 
traversing the grain. For small molecules the energy released generally results in dissociation, but 
for hydrocarbons this most likely results in only the loss of an H atom, with the carbon skeleton 
remaining intact. The sticking coefficient is therefore expected to be 

Se{Z >0)^{1- Pes) (1 - e-'^/'^) « 0.5 (1 - e-'^/^^) . (30) 

Figure 7 shows measured Sg for electron recombination with positive hydrocarbon ions, together 
with fitting function (30). The measured sticking coefficients generally agree with equation (30) to 
within about a factor of 2 except for the value measured for naphthalene CioHg", which is smaller 
than predicted by about a factor of 4. Why this should be smaller is unclear. 
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3.2. Ions 

Ions of interest (H"'" and C"*") have large ionization potentials, so we assume that they have 
a high probability of seizing an electron if they arrive at the surface of a grain, whether charged 
or neutral. While this assumption would be invalid for grains which are already highly positively 
charged, in practice the rate of arrival of positive ions at the surface of such a grain is in any 
case negligible. We compute the contribution of ions to the charging rate using equation (26) with 
Si = 1. 



4. Ambient Conditions 
4.1. Radiation Fields 

For some calculations, we adopt a blackbody spectrum for the radiation field, with color 

temperature Tc and dilution factor w, so that u^, = 4:TtwBu{Tc)/c. It is convenient to characterize 
the radiation intensity by G = u^Id/''^Wa.h^ where u^^^ is the energy density in the radiation field 
between 6eV and 13.6 eV and u'^^^ = 5.33 x 10~^^ergcm~^ is the Habing (1968) estimate of the 
starlight energy density in this range. ^ The radiation is cut off at 13.6 eV . 

For the diffuse ISM, we adopt the average interstellar radiation field (ISRF) spectrum in the 
solar neighborhood, as estimated by Mezger, Mathis, & Panagia (1982) and Mathis, Mezger, & 
Panagia (1983): 

fO ,hiy >13.6eV 

3.328 X lO-'^ergcm-s (/iz// eV )"^-^^^^ , 11.2eV < hu < 13.6eV 
8.463 X 10-13 erg cm-3(/ij,/eV)~^ , 9.26 eV < ^i^ < 11.2 eV ; (31) 
2.055 X 10-1^ erg cm-3 {hu/eVf-^^^^ , 5.04 eV < hu < 9.26 eV 
L {Attu/c) E-=i WiB^iTi) , hu < 5.04 eV 

the dilution factors Wi and blackbody temperatures Tj are given in Table 1. The total energy density 
in the ISRF of equation (31) is u = 8.64 x lO^^^ ergcm"^, with nj?^ = 6.07 x 10"-^^ ergcm^^ in the 
6-13.6 eV interval, or G = 1.13. 

The spectrum-averaged absorption efficiency factor is 

,^ > _ jQsbsUudu 

\Qa.hs) = r 1 , (,32j 

where QabsTra^ is the absorption cross section. In Figure 8 we display (Qabs) the ISRF and 
blackbody spectra with various values of Tc. 



uu}f^^ 



®For comparison, the interstellar radiation field estimated by Draine (1978) has u = 8.93 x 10 erg cm ^ between 
6 and 13.6 eV, or G = 1.68. 
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4.2. Scaling Law 

The photoelectric emission is dependent on the ambient conditions, which can be characterized 
by the shape of the radiation spectrum, the gas temperature T, and one additional parameter, 
depending on the ratio G/rig, which we take to be G\/T /rie- 

Unless otherwise noted, we will display results for a blackbody spectrum with Tc = 3 x 10^ K, 
cut off at 13.6 eV . In Figure 9, charge distributions are plotted for carbonaceous grains with various 
sizes in various environments. In Figure 10, the average electrostatic potential, (C/), is plotted as 
a function of grain size for two values of T, 100 K and 1000 K, and five values of G^/T /ue, ranging 
from 10^ to 10*^ K^/^ 

cm . We also provide results for conditions appropriate for the cold neutral 
medium (ISRF, T = 100 K, Ue = 0.03 cm'^, G^/T/ue = SSOcm^K^S) and the warm neutral 
medium (ISRF, T = 6000 K, Ue = 0.03 cm-^, Gy/T/ue = 2900cm3KV2). Note that for given 
Gs/r/ue, the computed potentials show hardly any dependence on T. 

For low values of G^T/ue, the grains can be negatively charged. This is especially so for grains 
with a < le ^ 10 A, for which the photoemission rate scales approximately with rather than a^. 
However, once a < 4A, the potential rapidly rises to positive values; this is a consequence of the 
rapid decrease in Sg for grains with Z < as a decreases in this range (eq. 28). The wiggles in the 
curves for high GVT/ue are due to saturation at the positive threshold charge Zjp. As the grain 
size increases, fz{Zip) — > 1, until the size at which Zjp increases by 1 is reached. Thus, although 
(Z) increases monotonically with a, it plateaus until Zjp is incremented, so that {U) decreases. 
Once Zip does increment, {U) does not abruptly increase, because the photoemission rate depends 
on hu — hupet and hu — hupdt ; thus charge is gradually concentrated in the next higher charge state 
as a increases. 

To see why G-\/T/ Ue is a good parameter to describe grain charging, we make the approxima- 
tion that the charge state on any grain remains constant, i.e. 

Je ~ Jpe "I" <-^ion ■ (^^) 

We neglect Jion by virtue of the large proton mass and consider a positively charged grain with 
?7 ~ IV. The Draine &: Sutin (1987) expression for J for a grain with static dielectric constant 
e(0) > 1 and < is 

which expands, for eU/kT :$> 1, to give 
We express 

Jpe (X G (Qabs) a^giU, a) , (36) 



1 + 



-2^ 



1/2 



(34) 
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where (Qabs) is the average value of Qabs over the range of absorbed photon energies. Here g{U, a) 
must be a decreasing function of U and g depends on a through the size dependence of the yield 
Y (see equations 12 and 13). Thus, 

g{U, a) He 

Note that this result applies only when G/ue is large enough (for given T and a) to keep the 
grains strongly positively charged (with ell > kT). For low values of G/ue, Jion plays a significant 
role in the grain charging. This introduces an additional dependence on the ionization fraction 
X = Ue/nii, since for the lowest x (i.e. x < 2 x 10~^), the ions are predominantly C"*", whereas 
H"*" dominates for higher x, and Jion depends on the ion mass. We assume H+ in our calculations 
for the CNM and WNM and C"*" in our other calculations. For the values of G\fT juf. and T 
considered here, the dependence on ion mix is only significant for the smallest grains, and even 
then extreme changes in the ion mix lead to only modest changes in the results. In Figure 11 we 
display (G\/T'/ne)o, the value of the charging parameter for which {Z) = 0, for carbonaceous and 
silicate grains, a blackbody radiation field with Tc = 3 x 10^ K, and various gas temperatures. As 
a decreases beyond 4A it becomes very difficult for the grains to charge negatively, because Sg 
becomes very small for Z <0; thus the curves plunge sharply there. 

Incidentally, the presence of (Qabs) in equation (37) explains the maxima in (U) at a ~ 100 A 
(see Figure 10), since Qabs peaks for 27ra ^ A, and the photoelectric emission is primarily from 
radiation with A ~ 1000 A. 



5. Photoelectric Heating Efficiencies 

The gas heating rate per grain due to photoelectric emission is given by 

rU(«) = E /^(^) Ke,M + r;d(«)] • (38) 

z 

The contribution from the photoemission of valence electrons is 

T'^,^,{a) = 7ra2 / duYQ,^,-^ / dEfE{E)E , (39) 

where £"01111 = when Z >0 and is given by equation (7) when Z < 0, and E^ax = hu—hupet+Eynin- 
The photoelectron energy distribution fE{E) = fE{E)/y2 (see eq. 10 and the discussion following 
it). When Z < 0, the contribution from photodetachment is given by 

rpd(a) = / du(jpAt{t^)^{hv - hi^pdt + -Emin) • (40) 
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We compute the total efficiency for conversion of absorbed radiation into gas fieating, 

T'^M - A'^ja) 

er{a) = — — ^ ■ (41) 

7ra^CUrad(Qabs) 

Here Agj.(a) is tfie rate of energy removal from the gas due to the accretion of charged particles 
onto the grain, and is given by 

(ShT\ ^^'^ 
7ra^A{Ti,^i)kT , (42) 
TTTTij / 
I, ' 

where the sum runs over electrons and ions and Tj and are defined below equation (26) . Expres- 
sions for A can be found in Draine Sz Sutin (1987). 

Figures 12 through 14 show results for the gas heating efficiency for the same conditions for 
which potentials were displayed in Figure 10. The heating efficiency generally decreases as grain 
size a increases. For low values of Gs/T/ue, this results primarily from the size-dependent yield 
factor yi (equation 13); the heating efficiency levels off at a 10^ A, where yi levels off. For 
higher G\fT juf.^ high grain potentials quench the photoemission and result in much lower heating 
efficiencies; in these cases the smallest grains (with a < 100 A), with somewhat lower potentials (sec 
Figure 10), dominate the photoelectric heating. The heating efficiencies drop dramatically when 
a < 4 A, as Se(.^ < 0) decreases rapidly and the grains charge positively. 



6. Net Photoelectric Heating Rate 

In this section we integrate the photoelectric heating efficiency over grain size distributions 
to find net heating rates. Cardelli, Clayton, h Mathis (1989) found that the extinction curve - 
and therefore the grain size distribution — varies depending on the environment through which the 
starlight passes, and that the variation can be roughly parameterized by Ry = Ay j Eb-v ^ the 
ratio of visual extinction to reddening. For the diffuse ISM, Ry k, 3.1; higher values are observed 
for dense clouds. 

Emission in the 3 to 60/im range, presumably generated by grains small enough to reach 
temperatures of 30 to 600 K or more upon the absorption of a single starlight photon (see, e.g. Draine 
&; Anderson 1985), imply a population of very small grains (with a < 50A). The non-detection of 
the lOjLtm silicate feature in emission from diffuse clouds (Mattila et al. 1996; Onaka et al. 1996) 
appears to rule out silicate grains as a major component of the a < 15 A population (Li &: Draine 
2000). Emission features at 3.3, 6.2, 7.7, 8.6, and 11.3;um (see Sellgren 1994 for a review) have 
been identified as C-H and C-C stretching and bending modes in polycyclic aromatic hydrocarbons 
(Leger & Puget 1984), suggesting that the carbonaceous grain population extends down into the 
molecular regime. The C abundance^ 6c in the very small grain population is still uncertain; 



By "abundance" , we mean the number of atoms of an element per interstellar H nucleus. 
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comparison of the observed diffuse galactic infrared and microwave emission with detailed model 
calculations for grains heated by galactic starlight imply 6c ~ 2—6 x 10~^ (Li &; Draine 2000; Draine 
& Lazarian 1998a,b; Draine k Li 2000). 

Weingartner & Draine (2000, hereafter WDOO) obtained size distributions with various values 
of 6c which reproduce the average observed extinction for lines of sight characterized by Ry = 3.1, 
4.0, and 5.5. Some of their distributions, designated "case A", were constructed so as to minimize 
the use of C and Si. For Ry = 4.0 and 5.5, "case B" distributions use the same amount of C and Si 
as those with Ry = 3.1. Li & Draine (2000) found that the diffuse galactic infrared emission is best 
reproduced with the WDOO size distribution with Ry = 3.1 and 6c = 6x 10~^. For this distribution, 
the 2175 A hump in the extinction curve is entirely due to the ultrasmall grain population. If this 
is generally the case throughout the ISM, then we would have 6c = 4 x 10^^ when Ry = 4.0 and 
6c = 3 X 10~^ when Ry = 5.5. We will usually display integrated heating rates computed for the 
size distributions with these (largest allowed) values of 6c. 

In Figure 15 we plot the net total gas heating per H nucleus and per Habing flux, 

itot _^ / ipe i dng,^^ ^^^^ 



c.s 



Gnn G da 



as a function of CVT/ue for our favored distributions from WDOO with T = 100 K and Tc = 
3 X 10^ K. The sum is over the carbonaceous and silicate populations. At a sing le GVf/ue, the 
heating rates obtained for the several distributions vary by a factor of ~ 2-4, showing the sensitivity 
of Ftot to the grain size distribution. Of course, the highest heating rates are for the distribution 
with Ry = 3.1 and 6c = 6 x 10~^, which contains the largest population of very small grains. 
For comparison, we also plot the Bakes & Tielens (1994) result. When G-v/T/ng > lO^K^/^cm^, 
their curve lies above our (Ry = 3.1, 6c = 6 X 10~^) curve, even though they assumed a smaller 
population of very small grains; this is because they adopted larger electron sticking coefficients 
than we do. In Figure 16 we show Ttot/Gnu for dust exposed to the ISRF; here we consider 
typical diffuse cloud dust with Ry = 3.1, for two possible values of 6c = and 6 x 10~^, and gas 
temperatures of 100 and 6000 K. Note that the net heating when T = 100 K does not differ greatly 
from that displayed in Figure 15 and that the grains have a net cooling effect for low values of 
GVr/ne when T = 6000 K. 

The photoelectric heating rate for the WDOO grain size distributions is fairly well reproduced 
by the following function: 

r,, . r~ T'^(a/-^ia = 10--5!1g„„ (44) 



g,S "mm 



where T is in K and OVT/ng is in K^/^cm'^. Values for the seven parameters in equation (44) 
are given in Table 2, along with the largest fractional error, err, for 10 < T < 10^ K and 10^ < 
cm'^. For most cases shown, we consider a blackbody spectrum with Tc = 
3 X 10^ K; we also consider the ISRF for Ry = 3.1. 
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Bakes & Tielens (1994) neglect heating from grains with a > 100 A. We find that this is 
usually a good approximation; large grains contribute significantly only for low G\fT jiif. and size 
distributions with few small grains. This is demonstrated in Table 2, where we give the fraction 
hs of the total heating due to grains with a < 100 A, for OVT/ne = 10^ 

cm and T = 100 K. 

For example, for Ry = 3.1, 6c = 0.0, and Tc = 3 x 10^ K, 25% of the heating is contributed by 
grains with a > 100 A. Bakes & Tielens also neglect the heating contributed by silicate grains. 
We find that silicates contribute < 25% of the total heating, with the maximum occuring for the 
distributions with the fewest very small grains. 

The rate of cooling due to charged particle collisions is significant, compared with the pho- 
toelectric heating rate, when T > 10'^ K. The following approximation is fairly accurate when 
10^ < r < 10^ K and 10^ < OVT/ue < 10^ kV2 cm^: 

Ag, = 10-28 erg cm^ s"^ nennT^^°+^'/''^ exp {D2 + D^x - D^x^) , (45) 

where x = \n.{G^/T / Ue) , T is in K, and G\fT jrif^ is in K^/^ cm^. In Table 3 wc give the values of 
Di and the maximum fractional error, err, for the above range of T and G\fT juf,. 

7. Heating and Cooling in H II Regions 

7.1. Stromgren Spheres 

Lyman continuum radiation from hot stars photoionizes the surrounding gas, resulting in an 
H II region. Studies of the heating in H II regions have found that photoelectric emission from 
dust can be important, compared with photoionization of H (Maciel &; Pottasch 1982; Oliveira 
h Maciel 1986; Maloney, Hollenbach, k, Townes 1992). In these studies, it was assumed that the 
photoelectric yield Y and absorption efficiency factor Qabs a^re independent of grain size and photon 
energy. Here, we apply our more detailed photoemission model. 

In H II regions, the radiation field includes photons with hv > 13.6 eV . To see the possible 
importance of these photons, we calculate integrated photoelectric heating and recombination cool- 
ing rates for blackbody spectra with no upper cutoff' energy, adopting Tc = 3.5 and 4.5 x 10^ K. Of 
course, at any given location in an H II region, there will be a break in the spectrum at 13.6 cV , due 
in part to absorptions along the path to the star and in part to the break in the stellar spectrum 
itself. Thus, spectra with no break and with a cutoff at 13.6 eV should bracket the range applicable 
in H II regions. 

We take ng/nn = 1 and T = 9000 K for the gas. In order to estimate likely G/nn values, 
we consider an 09 V exciting star and a point at the half-mass radius (i.e. at the distance from 
the star for which half of the H II region volume is enclosed). The radius of the H II region (the 
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"Stromgren radius") is found by balancing ionizations against recombinations: 

--(ill)''' 

where A^Ly is the rate at which the star produces ionizing photons and the case B recombination 
coefficient a « 2.6x lO-^^Tf"-^ cm^ where = T/IO"^ K (Osterbrock 1989). Since Rs oc n^'^ , 

1 /3 

G/nn oc n^' . Wc take A'^Ly = 3.63 x 10^^ s~-^, luminosity L = 4.4 x lO'^^ergs^^, and effective 
temperature T^^ = 3.59 x 10"^ K (Vacca, Garmany, & Shuh 1996). For nn = 0.1 (10'^) cm^'^, 
G/nH ^ 0.5 (10.) cm^. 

In Tables 4 and 5, we give Fpe/GraH and Ag^- /Gnu for G/uh = 0.1, 1, and 10 cm'^, both with 
and without the cutoff at 13.6 eV , for several of the WDOO grain size distributions. The heating 
rate typically increases by a factor of a few when the cutoff is removed, and the cooling rate 

is modestly affected, due to changes in the charging. When G/n-u = 1.0 cm^, grains have a net 
cooling (heating) effect for radiation with(out) a cutoff at 13.6 eV ; grains always have a net heating 
effect when G/riu = 10. cm^. When G/riu < 0.1 cm^ the grains are negatively charged, so that 
Tpe/Gnu ~ constant, and Ag^/Griu oc (G/nu)^^, roughly. 

The net steady-state heating rate due to photoionization is, in the notation of Spitzer (1978), 

^ anl mkT^ - {x2/h)kT) ^ 1.5 x 10-^^n^T,-°-8 (^^^^^^^^^g^^) ergcm^s-^ (47) 

where {tp)kTc is the mean kinetic energy per photoelectron, and {x2/<p2)kT is the mean kinetic 
energy per recombining electron for case B recombination; {ip) « 1.067(Tc/10^ K)^-^''^ for 16000 < 
Tc < 64000 K, and X2/<p2 ~ 0.67r4-°-^^ (Spitzer 1978). Comparing with the entries in Tables 4 
and 5, we find that the net heating from dust can be comparable to or exceed the photoionization 
heating when G/nn > 5cm^. 

In an actual H II region, the grain size distribution could be a function of distance from the 
exciting star, and could differ substantially from those found by WDOO. The very small grains 
might be destroyed in the harsh ionizing environment and the large grains might drift away from 
the star. A detailed study of the contribution of photoelectric heating in H II regions would have 
to account for these effects. 



7.2. The Warm Ionized Medium 

The warm ionized medium (WIM), also known as the diffuse ionized gas (DIG), extends to 
distances \z\ > Ikpc above the Galactic midplane (see Reynolds 1990, 1993 for reviews). It is not 
yet clear how the WIM is ionized and heated. The only known source of ionization with adequate 
power to maintain the WIM is Lyman continuum radiation from O stars, but it is not clear how 
the ionizing photons can reach locations with high \z\, since atomic gas is highly opaque to these 
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photons. Observed emission line intensity ratios can be approximately reproduced by models in 
which photoionization by a dilute radiation field dominates (Domgorgen & Mathis 1994). However, 
such models fail to account for observed variations in line intensity ratios with \z\ (see Reynolds, 
Haffner, & Tufte 1999 and references therein). 

Haffner, Reynolds, &: Tuftc (1999) found that the line intensity variations could be due to 
variations in the gas temperature; specifically, they suggested an increase from T « 7000 K at 
\z\ « 500 pc to T 10^ K at \z\ « 1500 pc. Reynolds et al. (1999) noted that the variations 
are more generally correlated with the electron density, suggesting a supplemental heat source 
that dominates photoionization heating at low densities. They found that a supplemental heat 
source T ^ 10~^^nH ergs~^ would result in a temperature profile Td^l) which could account for the 
observed line intensity ratio variations, and suggested photoelectric heating by dust as a possibility. 

In Figure 17, we plot the Td^l) profile inferred by Reynolds et al. (1999) from observed values 
of the [N II] /H a line intensity ratio. In this section we will investigate whether the combination 
of heating by photoionization and photoemission from dust can reproduce this profile. Reynolds et 
al. find that the electron density profile can be approximated by 

ne(|z|) = 0.125T4°-^V"°-^exp(-|z|/lkpc)cm-3 , (48) 

where T4 = T/10^ K and / is the WIM volume filling fraction. Reynolds et al. considered two cases 
for f{\z\): f{\z\) = 0.2 and f{\z\) = 0.1 exp(|z|/0.75 kpc) (Kulkarni k Heiles 1987). We calculate 
temperature profiles for four different cases (A-D), as summarized in Table 6, employing the above 
two prescriptions for the filling factor and two grain size distributions, for Ry = 3.1 and be = and 
6 X 10~^. We take the net heating rate due to photoionization Fpi = 1.5 x 10~'^^nlT^^'^ erg cm^ 
(appropriate for ~ 3.5 x 10^ K; equation 47) and cooling rate (excluding grain collisional cooling) 
Agas = 3.0 x 10~^^T4 '^ng ergcm^ s~^, which approximates the cooling function for ionized gas given 
in Figure 3.2 of Oster brock (1989), and find T satisfying 

Fpi + Fpe = Agas + Agr . (49) 

We plot the temperature profiles for the four cases A-D in Figure 17; we have adjusted the 
values of G so that the curves roughly lie on top of the Reynolds et al. curve. In Figure 18, we 
plot the associated values of the charging parameter CVT/ue and the net grain heating Ftot/^H 
as functions of \z\; note that Ftot/'^-H is not constant, since T and OVT/ue vary with \z\. For 
each of the cases A-D, the temperature profiles are too shallow compared with the Reynolds et 
al. curve. One might expect steeper curves to result if G is increased, although this might require 
reducing Fpi/rig to unreasonably small values. However, this is not usually the case, because the 
photoelectric heating efficiency begins to drop rapidly when G^T/ug > lO^K^/^cm^. We can 
obtain a reasonable match to the Reynolds et al. curve by altering case D somewhat, so that 
Fpi = 0.7 X lO'^'^nlT^^-^ ergcm^ s~^ and G = 0.27 (case E).^ Thus, we can reproduce the Reynolds 



This could arise, for example, for Tc w 26000 K. 
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et al. temperature profile by adopting a large population of very small grains and the Kulkarni k, 

Hcilcs (1987) formula for the WIM filling factor as a function of \z\. Other heating mechanisms, 
such as the dissipation of interstellar turbulence (Minter &: Spangler 1997), could also contribute, 
but do not appear to be required to account for the observed temperatures. 

8. Summary 

We have presented a model for photoelectric emission from dust. Here we first summarize 
how our model differs from the recent study by Bakes k, Tielens (1994) and the uncertainties that 
remain. 

We have re-evalulated the photoemission threshold energies for small grains, based in part on 
an empirical approach employing laboratory data for ionization potentials and electron affinities of 
PAH molecules. For negatively charged grains, we consider both photoemission from the valence 
band and the removal of attached electrons, and set the photoemission and photodetachment 
threshold energies equal to the minimum excitation energy for which an electron can effectively 
tunnel across the Coulomb barrier. Bakes & Tielens take the emission threshold equal to the 
electron affinity, which is an underestimate. Motivated by the PAH hypothesis. Bakes k. Tielens 
adopt a disk geometry for the smallest grains. The geometry affects the ionization potential through 
the capacitance, but this effect is small. We assume spherical grains for simplicity, and thereby 
retain consistency with the calculation of the charged particle accretion rates and Jion- 

We treat the variation of the photoelectric yield with grain charge and the distribution of 
photoelectron energies in a consistent manner (eq. 11). Whereas Bakes k Tielens adopt a delta 
function for the distribution of photoelectron energies, we adopt a parabolic form (eq. 10). We 
adopt a yield function for carbonaceous grains which, like that of Bakes k Tielens, approximately 
reproduces the measured ionization yield of coronene for grains with a 4A, but which agrees 
better with experimental determinations of the bulk graphite yield. Still, our bulk yields substan- 
tially exceed the experimental results; see the discussion following equation (16). In computing the 
enhancement factor for the yield of small particles, we evaluate the photon attenuation length la as 
a function of wavelength (eq. 14), rather than adopting la = 100 A, independent of wavelength, as 
Bakes k Tielens do. Recall, however, that the expression for the size-dependent yield enhancement 
(eq. 13) is based on a simple theoretical model; laboratory investigation of photoemission from 
microscopic particles is still in its infancy. 

Recent experimental results have been used to obtain new estimates for electron sticking effi- 
ciencies, as a function of grain size, shown in Figures 6 and 7. We have applied our photoemission 
and coUisional charging model to compute the efficiency with which grains convert absorbed ra- 
diation into gas heating via the photoelectric effect, as a function of grain size (§5, Figures 12 
- 14). Bakes k Tielens (1994) derived photoelectric heating rates for H I regions by integrating 
their heating efficiencies over an MRN size distribution with the lower cutoff size for the graphite 
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population extended to 3.5 A. We integrate the net gas heating rate over the size distributions of 
Weingartner & Draine (2000) (Figures 15 and 16) and we provide fitting functions for the result- 
ing total photoelectric heating and recombination cooling rates (§6). In a separate paper, we will 
examine the consequences of these revised rates for the thermal structure of the ISM. 

Finally, we have applied our model to estimate rates for photoelectric heating by dust grains 
in H II regions (§7). We find that photoelectric heating by dust can be important in H II regions 
when G/nu > 1 cm^. Photoelectric heating in the warm ionized medium (WIM) might explain the 
observed variations in emission line intensity ratios with distance from the Galactic midplane, if the 
abundance of ultrasmall grains is large enough. The required abundance is equal to that adopted 
by Li Sz Draine (2000) to account for the diffuse galactic infrared emission. 

We have investigated a range of grain sizes, gas temperatures, radiation field color temper- 
atures, and ratios G/rie of radiation intensity to electron density, and selected results have been 
presented. Interested readers can find a FORTRAN routine that implements the heating and cooling 
approximations of equations (44) and (45) on the World Wide Web at www.cita.utoronto.ca/~weingart. 
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Fig. 1. — The potential confining electrons in a grain with charge Ze. Shaded regions show 
occupied energy levels and W is the work function. Upper panel: Z > ~1. The photoemission 
threshold photon energy hup^t = IP, the ionization potential. Middle panel: Z = —1. The extra 
electron occupies the lowest unoccupied energy level (LUMO) of the neutral grain, lying a distance 
EA{Z = 0) below zero (EA is the electron affinity). The ionization potential IP{Z = —1) and 
photodetachment threshold energy /ifpdt both equal EA{Z = 0). The photoemission threshold 
energy hv^ct = IPv, the valence band ionization potential. IPy — IP ~ i^^bg) the energy band gap 
in bulk material; the equality is not exact due to quantum shifts in energy levels with grain size. 
Lower panel: Z < —\. When an electron acquires an energy -Emin above zero, it can tunnel out of 
the grain. Thus, /ifpdt = IP{Z) + E-aAn and /if pet = -f-fV(-Z^) + -Emin- When Z < —1, the maximum 
value of the confining potential is ^|z+i|e^/a; the parameter ^|z+i| is derived by Draine Sz Sutin 
(1987). 
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eq.(2) 
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W(graphite) 
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20 40 60 

Fig. 2. — First and second ionization potentials IP{0) and IP{1) for cyclic aromatic hydro- 
carbons and fullerenes. Solid curve is eq.(2). Curve labelled BT94 is the IP estimate of 
Bakes & Tielens (1994). Data: Lias et al. (1988) for C14H10 anthracene, C26H16 hexacene, 
C30H14 dibenz[bc,hl]-coronene, C30H16 pyranthene, C34H18 tetrabenz[a,cdj,lm]-perylene; Tobita et 
al. (1994) for CioHg naphthalene, C12H8 biphenylene, C12H8 acenaphthlylene, C12H10 biphenyl, 
C12H10 accnaphthcnc, C12H10 2-vinylnaphthalene, C13H10 fluorcnc, C14H10 diphenylacctylcne, 
C14H12 9,10-dihydrophcnanthrene, C14H9N acridinc, CieHio pyrcnc CiaHio fluoranthcne, C18H12 
tetracene, C20H12 pcrylcnc, C22H14 pentacene, C24H12 coronene; Chen et al. (1999) for C38H20 
benz[42]; Lichtenbcrgcr et al. (1991), Lichtenberger et al. (1992), and Steger et al. (1992) for Ceo, 
C70 fullerenes. 
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Fig. 3. — Electron affinity for neutral PAH molecules and Cgo- Solid curve is eq.(4). Curve 
labelled BT94 is estimate of Bakes (1992), used by Bakes & Tielens (1994). Data: Chen et 
al. (1999) for CeHg benzene, C24H12 coronene, C26H16 hexacene, C30H14 dibenz[bc,hl]-coronene, 
C30H16 pyranthene, C34H18 tetrabenz[a,cdj,lm]-perylene, C38H20 benz[42]; Shiedt k, Weinkauf 
(1997) for C14H10 anthracene; Ruoff et al. (1995) for CioHg naphthalene, C14H10 phenanthrene, 
CieHio pyrene, C18H12 tctraccnc, C18H12 bcnzanthraccnc, C18H12 chryscnc, C20H12 perylene, 
C20H12 benzo(a)pyrene, C22H14 pcntaccnc, C22H14 dibcnz(a,j)anthraccnc; Chen ct al. (1997) for 
C26H12 diindenochrysene, C28H14 debenzo[a,g]corannulene; and Wang, Ding, &; Wang (1999) for 
Ceo fullerene. 
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Fig. 4. — Photoionization yield for coronene, as measured by Verstraete et al. (1990) (solid) and 
as calculated using our model (dashed). 
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Fig. 5. — Photoelectric yield Y for neutral graphite and silicate grains as a function of incident 
photon energy hu, for several values of the grain size a, as indicated. 
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Fig. 6. — Sticking coefficients Se for electron attachment to neutral molecules, as a function of Nq, 
the number of atoms other than H (mostly C). Diamonds: acridine C13H9N, fluoranthene CigHio, 
pyrenc CieHio, tetracene C18H12, and perylene C20H12 (Tobita et al. 1992). Triangles: anthracene 
C14H10 (Tobita et al. 1992; Canosa et al. 1994; Moustefaoui et al. 1998), Ceo (Smith et al. 1993), 
and C70 (Spanel & Smith 1994). Filled symbols are for T = 300K. Open symbols for Ceo and C70 
are for T = 500 — 4500K. Open symbols for anthracene are for T = 48 — 170K (Moustefaoui et al. 
1998). Solid curve: adopted empirical fit given by eq. (28). Broken curve: Sg adopted by Tielens 
(1993) and Dartois k d'Hendecourt (1997). 




Fig. 7. — Sticking coefficients Sg for electron recombination with positive ions. Hexagons: T = 
300K experimental results of Abouelaziz et al. (1993) for C3HJ, CsH^, CeH^, C/H^, and CiqH^. 
Squares: T = BOOK results of Lefifaoui et al. (1997) for alkanes CH+, C2H+, CsHf , C4H^, C5H+ , 
C6H+3, C7H+ , and C8H+ . Triangle: T = 300K results of Rowe et al. (1995) for C7H+ and C14H+ 
phenanthrene. Diamonds: T = 300K results of Rebrion-Rowe et al. (1998) for C4H^ C4H^]^ CsH^ 
C6H4" CeH^ CsHy . Solid curve: adopted empirical fit given by eq. (30); shaded region is within 
factor of two. 
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Fig. 8. — Absorption efficiency factors for neutral carbonaceous and silicate grains, averaged over 
the interstellar radiation field (ISRF) and blackbody spectra with indicated color temperatures (cut 
off at 13.6 eV ). The kink at a = 50 A results from the Li & Draine (2000) prescription for blending 
PAH and graphite optical properties. 
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Fig. 9. — Charge distributions for carbonaceous grains with a = 3.5, 5, 7, 10, and 15 A and six 
sets of ambient conditions. For the upper four panels, we adopt a blackbody radiation spectrum 
(cut off at 13.6 eV) with Tc = 3 x lO'^K and T = 10^ K; the values of CVT/rie, in K^/^ cm"' are 
indicated. The bottom two panels are for the warm neutral medium (WNM) and cold neutral 
medium (CNM). 



-36- 




a (A) 

Fig. 10. — Average electrostatic potential ([/) for a blackbody radiation field with Tc = 3 x 10^ K 
(cut off at 13.6 eV). The top (bottom) panel is for silicate (carbonaceous) grains; curves labelled 
CNM (WNM) are for cold (warm) neutral media; solid (short-dashed) lines indicate T = 100 K 
(1000 K). For the T = lOOK and 1000 K curves, values of GyfTjue. are indicated in K^^ cm^. 
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Fig. 11. — The value of the charging parameter for which (Z) = 0, for carbonaceous and sihcate 
grains, a blackbody radiation field with Tc = 3 X 10^ K (cut off at 13.6 eV), and various gas 
temperatures, as labeled. 
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a (A) 

Fig. 12. — Gas heating efficiency for carbonaceous grains for two gas temperatures: 100 K (solid) 
and 10^ K (sliort-dashed) , four values of G\fT jn,,., as labelled, and = 3 x 10^ K. (The radiation 
is cut off at 13.6eV.) Curves labelled CNM (WNM) are for cold (warm) neutral media and the 
ISRF. 




Fig. 13. — Same as Figure 12, but for silicate grains. The curve for T = 1000 K and OVT/ng = 
10^ K^/^ cm^ lies entirely below the plot region. 
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Fig. 14. — Gas heating efficiency for carbonaceous grains, T = 100 K, GVr/ne = IO^kVs 
and four values of as indicated. (The radiation is cut off at 13.6 eV .) 
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Fig. 15. — The photoelectric heating rate minus the coohng rate due to gas phase ions and electrons 
colliding with and sticking to the grains; for a blackbody radiation spectrum (cut off at 13.6 eV ) with 
Tc = 3 X 10^ K, gas temperature T = 100 K, and various grain size distributions from Weingartner 
& Draine (2000), as indicated. For comparison, we also display the Bakes & Tielens (1994) result 
(B&T). 
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Fig. 16. — Same as Figure 15, but for the ISRF of equation (31), Ry = 3.1, and various values of 
be and gas temperature T. 



-43- 



12000 I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




6000 ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' 

0.5 1 1.5 2 

|z| (kpc) 

Fig. 17. — WIM temperature profiles as a function of distance from the Galactic midplane \z\. The 
solid curve was inferred by Reynolds et al. (1999) from observations of the [N II]/H a line intensity 
ratio; the other curves are results of calculations including photoelectric heating, for various sets of 
parameter values as indicated in Table 6. 
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Fig. 18. — Upper panel: Values of the charging parameter as a function of \z\ for the five cases 
of Table 6. Lower panel: The ratio of the net grain heating rate (i.e. photoelectric heating minus 
collisional cooling) to the gas density as a function of \z\, for the five cases of Table 6. 
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Table 1. Blackbody Components of the ISRF (Mathis et al. 1983) 
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Wi 


Ti/K 


1 


10-14 


7500 
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1.65 X 10-13 


4000 


3 


4 X 10-12 


3000 
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Table 2. Photoelectric Heating Parameters^ 





be 


case 


rad field*^ 


Co 






Co 




Ca 




Cfi 


err 




3.1 


0.0 


A 


BO 


5.56 


1.82 X 10- 


-2 


0.00492 


0.03368 


0.557 


0.666 


0.532 


0.15 


0.75 


3.1 


2.0 


A 


BO 


3.41 


4.27 




0.03700 


0.00569 


0.102 


0.494 


0.669 


0.19 


0.89 


3.1 


4.0 


A 


BO 


10.59 


4.05 X 10" 


-3 


0.01234 


0.01391 


0.808 


0.580 


0.573 


0.19 


0.94 


3.1 


6.0 


A 


BO 


3.66 


7.66 




0.00661 


0.01552 


0.094 


0.697 


0.482 


0.20 


0.96 


4.0 


0.0 


A 


BO 


4.30 


1.79 X 10- 


-4 


0.00572 


0.02488 


1.026 


0.701 


0.505 


0.15 


0.77 


4.0 


2.0 


A 


BO 


5.03 


2.27 




0.07441 


0.00345 


0.140 


0.412 


0.737 


0.18 


0.92 


4.0 


4.0 


A 


BO 


3.28 


4.44 




0.00786 


0.01219 


0.102 


0.686 


0.500 


0.20 


0.96 


5.5 


0.0 


A 


BO 


3.67 


2.42 X 10- 


-3 


0.07034 


0.00360 


0.725 


0.365 


0.814 


0.15 


0.77 


5.5 


1.0 


A 


BO 


5.45 


2.44 X 10- 


-2 


0.09594 


0.00225 


0.534 


0.418 


0.758 


0.18 


0.89 


5.5 


2.0 


A 


BO 


4.87 


7.28 X 10' 


-3 


0.00812 


0.01969 


0.656 


0.620 


0.534 


0.19 


0.94 


5.5 


3.0 


A 


BO 


5.93 


9.26 X 10- 


~2 


0.00675 


0.01648 


0.417 


0.683 


0.489 


0.20 


0.96 


4.0 


0.0 


B 


BO 


4.85 


6.59 X 10- 


-2 


0.00437 


0.04234 


0.415 


0.664 


0.508 


0.16 


0.86 


4.0 


2.0 


B 


BO 


5.97 


0.287 




0.00736 


0.02065 


0.302 


0.637 


0.528 


0.18 


0.93 


4.0 


4.0 


B 


BO 


7.97 


0.175 




0.00277 


0.03820 


0.382 


0.771 


0.401 


0.19 


0.96 


5.5 


0.0 


B 


BO 


4.04 


0.653 




0.11310 


0.00259 


0.198 


0.348 


0.809 


0.18 


0.89 


5.5 


1.0 


B 


BO 


3.38 


0.608 




0.00396 


0.02568 


0.192 


0.752 


0.439 


0.18 


0.93 


5.5 


2.0 


B 


BO 


1.10 


3.30 




0.00941 


0.01431 


0.087 


0.643 


0.524 


0.20 


0.95 


5.5 


3.0 


B 


BO 


3.47 


1.77 




0.00494 


0.01953 


0.140 


0.723 


0.456 


0.20 


0.97 


3.1 


0.0 


A 


ISRF 


4.35 


2.63 X 10- 


-5 


0.00242 


0.03003 


1.235 


0.827 


0.399 


0.15 


0.73 


3.1 


2.0 


A 


ISRF 


4.87 


0.948 




0.02576 


0.00766 


0.188 


0.490 


0.650 


0.18 


0.88 


3.1 


4.0 


A 


ISRF 


7.41 


0.772 




0.03895 


0.00606 


0.239 


0.408 


0.708 


0.19 


0.94 


3.1 


6.0 


A 


ISRF 


9.30 


0.248 




0.00697 


0.01848 


0.365 


0.633 


0.509 


0.20 


0.96 



^See eq (44). 

^BO refers to a blackbody spectrum with Tc = 3 x 10^ K and the ISRF of Mathis et al. (1983) is 
defined in eq (31). 

'^Fraction of the total heating due to grains with a < 100 A, for G\/T/ne = lO^K^/^cm^ and 
T= 100 K. 
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Table 3. Collisional Cooling Parameters^ 



Rv 


be 


case 


rad field*^ 


Do 




D2 


^3 


L>4 


err 


3.1 


0.0 


A 


BO 


0.4800 


1.783 


-7.617 


1.655 


0.06326 


AIT 

0.17 


3.1 


2.0 


A 

A 


BO 


0.5322 


1.325 


-6.047 


1.435 


0.05378 


0.14 


3.1 


4.0 


A 

A 


BO 


0.4336 


2.108 


-7.359 


1.736 


0.06502 


0.14 


3.1 


6.0 


A 

A 


BO 


0.4270 


2.120 


-7.301 


1.786 


0.06731 


0.15 


4.0 


0.0 


A 


BO 


0.5332 


1.379 


-7.125 


1.504 


0.05732 


0.17 


4.0 


2.0 


A 

A 


BO 


0.4425 


2.019 


-7.755 


1.712 


f\ f\/^ AC\ A 

0.06424 


0.14 


4.0 


4.0 


A 

A 


BO 


C\ ACi Ai'7 

0.4847 


1.493 


-6.109 


1.528 


0.05751 


A 1 O 

0.13 


5.5 


0.0 


A 


BU 


0.4184 


2.129 


-9.082 


1.818 


A A£?A A O 

0.06943 


A 1 

O.lo 


5.5 


1.0 


A 


BO 


0.4981 


1.510 


-7.237 


1.535 


0.05777 


0.14 


5.5 


2.0 


A 


BO 


0.4851 


1.501 


-6.775 


1.533 


0.05762 


0.13 


5.5 


3.0 


A 


BO 


0.4948 


1.442 


-6.377 


1.505 


0.05645 


0.13 


4.0 


0.0 


B 


BO 


0.4568 


1.876 


-7.810 


1.688 


0.06412 


0.15 


4.0 


2.0 


B 


BO 


0.4597 


1.594 


-6.675 


1.599 


0.06094 


0.14 


4.0 


4.0 


B 


BO 


0.5155 


1.340 


-5.791 


1.449 


0.05411 


0.13 


5.5 


0.0 


B 


BO 


0.4360 


1.996 


-8.372 


1.735 


0.06572 


0.15 


5.5 


1.0 


B 


BO 


0.5140 


1.349 


-6.676 


1.472 


0.05536 


0.14 


5.5 


2.0 


B 


BO 


0.4633 


1.690 


-7.158 


1.607 


0.06055 


0.14 


5.5 


3.0 


B 


BO 


0.4677 


1.780 


-7.173 


1.615 


0.06027 


0.13 


3.1 


0.0 


A 


ISRF 


0.4291 


2.406 


-8.357 


1.714 


0.06354 


0.15 


3.1 


2.0 


A 


ISRF 


0.5232 


1.678 


-5.942 


1.339 


0.04813 


0.14 


3.1 


4.0 


A 


ISRF 


0.3959 


2.380 


-6.554 


1.575 


0.05674 


0.13 


3.1 


6.0 


A 


ISRF 


0.3632 


2.937 


-7.601 


1.742 


0.06228 


0.12 



^See eq (45). 

"^BO refers to a blackbody spectrum with Tc = 3 x lO^K and the ISRF of Mathis et al. (1983) is 
defined in eq (31). 
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Table 4. Photoelectric Heating Rates for H II Regions 





lO^bc'' 


case 


Tc^ G/nH = 


0.1 s 


1.0 


10. 


0.1 


1.0 


10. 


3.1 


0.0 


A 


3.5 


1.03 


0.83 


0.40 


3.58 


2.91 


1.25 


3.1 


6.0 


A 


3.5 


2.35 


2.06 


1.16 


8.15 


7.07 


3.59 


4.0 


4.0 


A 


3.5 


1.54 


1.36 


0.76 


5.37 


4.66 


2.35 


5.5 


3.0 


A 


3.5 


1.09 


0.97 


0.55 


3.83 


3.34 


1.68 


4.0 


4.0 


B 


3.5 


1.57 


1.38 


0.77 


5.49 


4.76 


2.38 


5.5 


3.0 


B 


3.5 


1.06 


0.94 


0.54 


3.74 


3.28 


1.64 


3.1 


0.0 


A 


4.5 


1.13 


0.91 


0.43 


5.93 


4.73 


1.80 


3.1 


6.0 


A 


4.5 


2.57 


2.26 


1.24 


13.1 


11.1 


5.02 


4.0 


4.0 


A 


4.5 


1.68 


1.48 


0.82 


8.62 


7.31 


3.25 


5.5 


3.0 


A 


4.5 


1.19 


1.06 


0.59 


6.13 


5.21 


2.31 


4.0 


4.0 


B 


4.5 


1.72 


1.51 


0.83 


8.81 


7.46 


3.30 


5.5 


3.0 


B 


4.5 


1.16 


1.03 


0.58 


5.96 


5.08 


2.26 



^10 ^^ergs ^ 

''Blackbody spectrum cut off at 13.6 eV 
'^Pull blackbody spectrum 

'^Rv = Ay/Eb-v, ratio of visual extinction to reddening 
•^C abundance in very small grain population 

flO^K 

^ Value for G/nji, in cm^ 
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Table 5. Recombination Cooling Rates for H II Regions 





be' 


case 


Tc^ G/nn = 




A/Gn-a 


a,b 




A/Gn-E 


a,c 


0.1 ^ 


'• 1.0 


10. 


0.1 


1.0 


10. 


3.1 


0.0 


A 


3.5 


7.43 


1.02 


0.22 


7.91 


1.36 


0.36 


3.1 


4.0 


A 


3.5 


29.4 


3.47 


0.64 


30.2 


4.17 


1.00 


4.0 


2.0 


A 


3.5 


18.7 


2.20 


0.41 


19.2 


2.66 


0.64 


5.5 


1.0 


A 


3.5 


13.2 


1.55 


0.29 


13.6 


1.87 


0.45 


4.0 


2.0 


B 


3.5 


18.9 


2.23 


0.42 


19.4 


2.69 


0.65 


5.5 


1.0 


B 


3.5 


13.1 


1.53 


0.28 


13.5 


1.84 


0.44 


3.1 


0.0 


A 


4.5 


7.46 


1.05 


0.23 


8.29 


1.58 


0.44 


3.1 


4.0 


A 


4.5 


29.4 


3.51 


0.67 


30.9 


4.68 


1.21 


4.0 


2.0 


A 


4.5 


18.7 


2.23 


0.43 


19.6 


2.98 


0.78 


5.5 


1.0 


A 


4.5 


13.2 


1.57 


0.30 


13.9 


2.10 


0.55 


4.0 


2.0 


B 


4.5 


18.9 


2.26 


0.43 


19.9 


3.02 


0.79 


5.5 


1.0 


B 


4.5 


13.1 


1.55 


0.29 


13.8 


2.06 


0.54 



'^lO"^^ ergs~^ 

''Blackbody spectrum cut off at 13.6 eV 
'^FuU blackbody spectrum 

'^Rv = Ay/Eb-v, ratio of visual extinction to reddening 
abundance in very small grain population 

^lO^K 

s Value for G/nn, in cm^ 
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Table 6. WIM Heating Examples 



Case 


Grain Size Distribution 


WIM Filling Factor 


rp./n^T,-°-« 


G 








10 erg cm^ s ^ 




A 


Rv = 3.1, he = 0.0 


0.2 


1.5 


0.3 


B 


Rv = 3.1, he = 0.0 


0.1exp|z|/750pc 


1.5 


0.3 


C 


Rv = 3.1, 6c = 6 X 10-5 


0.2 


1.5 


0.14 


D 


Rv = 3.1, 6c = 6 X 10-5 


0.1exp|z|/750pc 


1.5 


0.1 


E 


R\ = 3.1. he = 6 X 10-'* 


0.1 ~ /750 pc 


0.7 


0.27 



